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Abstract 

Open-flavor strong decays are mediated by qq pair production, which 
is known to occur dominantly with '^Pq quantum numbers. The relation of 
the phenomenological ^Fq model of these decays to "microscopic" QCD decay 
mechanisms has never been clearly established. In this paper we investigate qq 
meson decay amplitudes assuming pair production from the scalar confining 
interaction (sKs) and from one gluon exchange (OGE). sKs pair production 
predicts decay amplitudes of approximately the correct magnitude and D/S 
amplitude ratios in bi lott and oi pir which are close to experiment. The 
OGE decay amplitude is found to be subdominant in most cases, a notable 
exception being '^Pq — > ^Sq + ^Sq. The full sKs + OGE amplitudes differ 
significantly from ^Pq model predictions in some channels and can be distin- 
guished experimentally, for example through an accurate comparison of the 
D/S amplitude ratios in bi — > ujtt and ai pn. 
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I. INTRODUCTION: THE ^Pq MODEL. 

A. Summary of the ^Pq model. 

Strong decays constitute a rather poorly understood area of hadron physics. A phe- 
nomenological model of hadron decays was developed in the 1970s by LeYaouanc et al [|T]-§| , 
which assumes as earlier suggested by Micu 0] that during a hadron decay a qq pair is pro- 
duced from the vacuum with vacuum quantum numbers, J^^ = 0^^. Since this corresponds 
to a '^Pq qq state, this is now generally referred to as the ^Pq decay model. The ^Pq pair 
production Hamiltonian for the decay of a qq meson A to mesons B + C is usually written 
in a rather complicated form with explicit wavefunctions ||5|-§, which in the conventions of 
Geiger and Swanson ^ (to within an irrelevant overall phase) is 

{BCmA) =^ JJ e^^-- vl/^(f) {a),, ■ (^V^ + ^Vc + Pb) + y) ^ci^ - y) 

S{Pa-Pb-Pc) (1) 

for all quark and antiquark masses equal. The strength 7 of the decay interaction is regarded 
free constant and is fitted to data [H. 
Studies of hadron decays using this model have been concerned almost exclusively with 
numerical predictions, and have not led to any fundamental modifications. Recent studies 
have considered changes in the spatial dependence of the pair production amplitude as a 
function of quark coordinates but the fundamental decay mechanism is usually not 
addressed; this is widely believed to be a nonperturbative process, involving "flux tube 
breaking" . There have been some studies of the decay mechanism which consider an alter- 
native phenomenological model in which the qq pair is produced with ^Si quantum numbers 
p!0| ; this possibility however is rejected by experiment 0. In this paper we first review the 
phenomenological ^Pq model and its predictions and develop a diagrammatic representa- 
tion, following which we will test candidate QCD (quark-gluon) decay mechanisms for these 
strong decays. 
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An equivalent formulation of the ^Pq model regards the decays as due to an interaction 
Hamiltonian involving Dirac quark fields, 

Hi = g J d^x , (2) 

which in the nonrelativistic limit gives matrix elements identical to (1) with the identification 

7=^ (3) 
' 2mq ^ ^ 

where ruq is the quark mass of the produced pair. The operator gipip leads to the decay 

{qq)a {<1(1)b + {(l<l)c through the b^S term. 

This model makes no reference to color, which if included would simply change the 
definition of the interaction strength 7; since 7 is fitted to data, this would not change the 
predictions for meson decays. An ad hoc feature of the ^Pq model is to allow only diagrams 
in which the qq pair separate into different final hadrons. This was originally motivated by 
experiment (specifically the weakness of — > pvr); in the QCD-based current-current decay 
models that we shall discuss subsequently the absence of these hairpin diagrams is a natural 
consequence of the production of the qq pair in a color octet state. 

To determine a decay rate we evaluate the matrix element of the decay Hamiltonian, 
which is of the form 

{BC\Hi\A) = hf, S{A - B - C) . (4) 

This hfi decay amplitude can be combined with relativistic phase space to give the 
differential decay rate, which is 

dTA^BC o PEbEc |2 

= ZTT \hfi\ (5) 

where we have set A = and P = \B\ = \C\. An equivalent result is quoted by Geiger and 
Swanson 0] as the "actual phase space" case of their equation (20). We shall usually quote 
results for the amplitude MlbcSbc^ defined by 

rA^Bc = 2n^^Y.\MLs\' . (6) 

Ma ls 



For 5* = this M.ls is the coefficient of Ylm in /i/j, and in general is the amphtude 
{JM{Lbc, Sbc)\BC) in the BC final state. 

We have developed a diagrammatic technique for evaluating the hfi matrix element, 
which is discussed in detail for p —>■ nn in the ^Pq model in Appendix A. For qq meson 
decays there are two independent diagrams, shown below. 




Figure 1. The two independent qq meson decay diagrams in the ^Pq decay model. 



B. ^Pq model decay rates for light hqq = and hqq = 1 mesons. 

Since the ^Pq model involves a free phenomenological interaction strength 7, comparison 
with experiment requires a fit to several rates (to determine both 7 and the wavefunction 
parameters) or a determination of amplitude ratios in decays with more than one partial 
wave, in which 7 cancels. Here we illustrate both applications; first we will evaluate the 
dominant two-body decay rates of light nonstrange mesons with Lgg = and Lgq = 1, and 
subsequently we will evaluate D/S ratios in the decays 61 un, ai pvr and hi —>■ pvr. 

An M.LS decay amplitude in the ^Pq model with SHO wavefunctions is proportional to 
a polynomial Vls{x) in x = P/P times an exponential, 

-M.5 = ^^T7^P.5(x)e--^/-. (7) 

For the cases considered here these polynomials are 
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P^"-'"^'"'--^-' (10) 



p,3p,-.s,+.s., = +^ (1 _ 2,.) (13) 



- °" = -§575 (1 - r J (") 

p,;p.-.s.+.s., = _|!_ . (15) 

For physical decays there is an additional multiplicative flavor factor, as discussed in 
Appendix A. In figure 2 we show these decay rates for physical mesons for a wide range of 
wavefunction length scales /3 (defined in App.A). 




Figure 2. Representative light qq meson decay rates in the ^Pq decay model. Parameters j3 
variable, 7 = 0.5. 
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We set 7 = 0.5 in the figure because this is near the optimum fitted value. A value of 
P ~ 0.35-0.4 GeV evidently gives a reasonably accurate description of the relative decay 
rates (see also Table 1 below). Many recent quark model studies of meson and baryon 
pJ| decays use a value of (3 = 0.4 GeV, at the upper end of this range. 

A simple estimate of the optimum f3 for these SHO wavefunctions follows if we require 
maximum overlap with Coulomb plus linear wavefunctions \E'. With typical quark model 
parameters of = 0.6, b = 0.18 GeV^ and niq = 0.33 GeV, the overlap | / '^*ipsHO (i^xp 
is maximum at (3 = 0.316 GeV for Lqg = (99.4% overlap) and (3 = 0.274 GeV for 
Lgq = 1 (99.6% overlap). Curiously, this /5 would imply a narrow /o(1300) in the ^Pq 
decay model. Comparison with figure 2 shows that these f3 estimates are similar to the val- 
ues required to describe decays, although there is some evidence of a systematic discrepancy, 

P decay / Pcou.+lin. ~ 1.2. 



Table 1. ^Pq Model Fit to Light Meson Decay Rates. 


Decay 


Expt. [0 


■^Pq Theory 


D/S Expt |12,13| 


D/S =^Po Theory 


P — i> TTTT 


151. MeV 


79. MeV 








157. MeV 


170. MeV 






02 — >• piT 


72. MeV 


54. MeV 






ai pn 


400. MeV 


545. MeV 


-0.09(2) 


-0.154 


hi un 


142. MeV 


143. MeV 


+0.260(35) 


+0.292 


hi — *• pn 


360. MeV 


383. MeV 




-0.111 


K*(1430) ^ Ktt 


287. MeV 


166. MeV 








150-400 MeV (not fitted) 


271. MeV 







In Table 1 we show the result of a fit to a representative set of seven well established light 
S- and P-wave qq meson decay rates, with j3 and 7 taken as free parameters. The masses 
assumed were = 0.138 GeV, Mr = 0.496 GeV, Mp = 0.77 GeV, = 0.782 GeV, 
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Mh^ = 1.17 GeV, = 1.23 GeV, M^, = 1.231 GeV, M/^ = 1.3 GeV, M„, = 1.318 GeV 
and Afx* = 1.429 GeV. All but Kq(1430) — > Ktt are nonstrange systems, which we chose 
to avoid parameter differences due to strange quarks. The Kq(1430) is included because we 
found that it is quite important to include a light ^Pq qq decay, since the rate ^Pq — >^So+^So 
is quite sensitive to (3, and the Kq(1430) is the only well established ^Pq resonance. On 
fitting these widths using Eqs.(6-15) (minimizing Yld=iiX*^AtBC I^'^a^bc ~ 1)^)' the 
parameters 

(3 = 0.397 GeV (16) 

7 = 0.506 . (17) 

Clearly the most important discrepancy in the table is p tttt, which is well known to be 
a problem relative to the decays of P-wave qq mesons in this model. The Kg (1430) width is 
also rather smaller than experiment. 

If we attempt constrained fits with fixed (3, we find serious disagreement with experiment 
for even moderate changes away from (3 — 0.4 GeV. For example, Bi (3 — 0.3 GeV we find 
an /2 broader than the ai and hi and a Kq(1430) width near zero. Conversely, if we increase 
l3 to 0.5 GeV we find quite small partial widths for p —>■ tttt, a2 — * pvr and /2 tttt relative 
to the bi, Oi and hi, on average only 0.31 of bi — > cutt; experimentally the ratio is 0.89. 

The interesting question of the width of an /o(1300) qq state is unfortunately problemat- 
ical in the ^Pq model, since like Kq(1430) — ^ Ktt this rate has a node near P = 0.30 GeV, and 
increases rapidly with /? above the node. For our fitted parameters we find ^ fo(i3oo)-,-jvn- — 
271 MeV, but comparison of the predicted and observed Ko(1430) widths suggests that the 
actual /o(1300) partial width to tttt is ~ 450 McV. We shall subsequently show that the 
transition ^Pq ^^So+^Sq has an unusually large OGE decay amplitude, which may explain 
why the hght scalars /o(1300) and Kq(1430) are broad states despite the node in the ^Pq 
decay amplitude. 
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C. ^Pq results for D/S ratios 



Sensitive tests of decay models are possible in decays with more than one partial wave, 
because one can measure the relative phases as well as the magnitudes of the decay ampli- 
tudes in these processes. 

Multiamplitude decays require at least one final meson to have nonzero spin. For the 
decays considered here the relevant final states are pvr and uott; other possibilities are excluded 
by phase space. These final states have = 1'*', so of all the light L^g = and L^g = 1 
mesons only the l"*" mesons have more than one partial wave; these decay to both S- and D- 
wave vector + pseudoscalar final states. On increasing L^^ we next encounter multiamplitude 
decays in the D states, for example in — > pvr (P,F) and Ti2 /27r (S,D,G). 

Here we will calculate the D/S amplitude ratios in the decays bi uotx and ai pvr. 
D/S for /ii ^ pTT, which has not been measured, is theoretically equal to hi ujtt to within 
small phase space differences, since these are both ^Pi ^^Si+^Sq decays. The bi is clearly 
the most attractive of these experimentally, due to the smaller widths of the bi and u. 

We assume an initial polarization ^^^(J^ = +1) and determine the hfi matrix element to 
a final state with a specific uj polarization that allows both L values; we choose uj{Sz = +1). 
Using the diagrammatic techniques of Appendix A, we find that the two diagrams di and 
d2 give equal contributions, and the total hfi is 

7 2^ f/, 2 , 2 



h 



fi 



:{(l - '-x^ ) Y,,in) + ^x^y.o(l^)} e-^^/- . (18) 



,+ (+,)^^(+,)^+ 7rV4/3i/2 35/2 V 

This matrix element determines the relative S and D amplitudes in the \uj7c) final state; 
since the = 1+ |u;7r) state is of the form 

luon) = as\'%) + aofBi) (19) 

it follows from a Clebsch-Gordon decomposition that the amplitude to find a (+i)-polarized 
u in an un pair with decay direction Q must be 



We can therefore read the ao/cLs ratio directly from the hfi matrix element (18), which gives 



as 



X 



bl — ^iUTT 



23/2 



32 (1-1x2) 



(21) 



This is equivalent to the result (4.11,12) of LeYaouanc et al. [Q] when one specializes their 
result to Rh^ = R^j = Rn, which is our l/j3. This D/S ratio is also implicit in (14,15). Since 
the strength of the pair production amplitude and some of the momentum dependence in 
the overlap integrals cancel out in this ratio, there is less systematic uncertainty than in the 
decay rates. Proceeding similarly for ai pvr we find 



as 



21/2 



X 



ai^pn 32 (1 g X^) 

The ratio of D/S ratios in these decays is especially interesting in the ^Pq model because all 
the dependence on the spatial wavefunctions cancels. We then find 



ap 




a-s 


ai— >p7r 


ap 




as 





(23) 



Experimentally these ratios are 



ap 
as 



+0.260 ±0.035 



bl — >uJ7r 



(24) 



and from a recent ARGUS measurement ||13| 



so the ratio is 



ap 
as 



a-p 



-0.09 ±0.02 



(25) 



ai— *p7r 



ap 




as 





-0.35 ±0.09 . 



(26) 



Theory and experiment for these D/S ratios are compared in figure 3; a best fit to these D/S 
ratios alone gives (3 = 0.448 GeV, for which (D/S)^,^^^ = ±0.219 and (D/S)^,^^^ = -0.115. 
(The small theoretical departure from —1/2 is due to phase space differences.) 
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The experimental ratio (26) is approximately consistent with the ^Pq prediction but 
shows the need for a more accurate experimental determination, especially in ai —>■ pn. 
This important D/S ratio has only been measured in one experiment; an improved value 
would allow an interesting test of the ^Pq model, since the OGE decay mechanism predicts 
a departure of the ratio (23) from —1/2; we will discuss this in the next section. 

In summary, the decay rates of light Lqq = and Lqq = 1 mesons (figure 2) support 
j3 ~ 0.40 GeV, whereas the D/S amplitude ratios in bi un and ai — » pn (figure 3) suggest 
a larger value, f3 ^ 0.45 GeV. These are both significantly larger than the (3 ~ 0.3 GeV 
which gives the maximum overlap with Coulomb plus linear wavef unctions. 
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Figure 3. D/S amplitude ratios in bi — > un and ai pn with ^Pq model predictions. 



II. "MICROSCOPIC" DECAY MODELS 
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A. Candidate pair production Hamiltonians. 



In microscopic decay models one attempts to describe hadron strong decays in terms of 
quark and gluon degrees of freedom. The quark-gluon decay mechanism should give similar 
predictions to the reasonably accurate ^Pq model, and should determine the strength 7 of the 
^Po interaction in terms of fundamental QCD parameters. There has been little previous 
work in this area. One exception is the study of open charm decays of cc resonances by 
Eichten et al. [|I^, who assumed that decays are due to pair production from the confining 
interaction. We will confirm this assumption for light quarkonia in most but not all cases. 
Eichten et al. however assumed that confinement is a Lorentz vector interaction; this is 
now believed to be incorrect, and the currently accepted scalar form leads to quite different 
matrix elements. 

We begin by assuming that these strong decays are driven by the same interquark Hamil- 
tonian which determines the spectrum, and that it incorporates scalar confinement and one 
gluon exchange. These interactions and their associated decay amplitudes are undoubt- 
edly all present and should be added coherently. We will determine decay rates for the 
two-body decays discussed in the previous section, assuming the appropriate field-theoretic 
generalization of the usual interquark Hamiltonian. 

The current-current interactions due to the scalar confining interaction and OGE can be 
written in the generic form 



The current J" in (27) is assumed to be a color octet. The currents J (with the color 
dependence A"/2 factored out) and the kernels K{r) for the interactions are 




(27) 



scalar confining interaction 



J = < 



color Coulomb OGE 



(28) 



('?/'7'?/')t transverse OGE 
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K 



+ jlbr + So] scalar confining interaction 



color Coulomb OGE 
transverse OGE . 



(29) 



+as/r 
-Ois/r 

The scalar confinement kernel is normalized so that hr + 5*0 is the static scalar potential of 
a color-singlet qq pair. 

We refer to this general type of interaction as a JKJ decay model, and to the specific 
cases considered here as sKs, j'^Kj° and j^Kj"^ interactions. The decay Hamiltonian assumed 
by Eichten et al. in their Eq. (3.2) is a special case of our form (27) with a j^Kj^ 
interaction. 



B. JKJ A^BC matrix elements: general results. 

The A— >^BC decay matrix element of the JKJ Hamiltonian (27) involves a pair-production 
current matrix element {qq\J\^) times a scattering matrix element {qf\J\qi)- Diagrammati- 
cally this corresponds to an interaction between an initial line and the produced pair, and 
there are four such diagrams (see figure 4). We label these using the ^Pq quark line diagram 
labels di and (i2, with an additional g or g subscript denoting which initial line the produced 
pair interacts with. 

We specialize to the diagram dig for illustration. As in the ^Pq model (Appendix A) 
there is a Fermi "signature" phase, a flavor factor, and a spin+space overlap integral, 

(-SC* I -ff/ 1 j4^(^-i^^ -^signature ' -^color ' -^flavor ' Ispin+space (^^) 

where 

^spin+space ^spin+space '^(^ B C ) • (^-^) 
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iq 



2q 





".5 



Figure 4. The four independent qq meson decay diagrams in JKJ decay models. 

The signature is again /signature = {—^T because there are three hne crossings. A new- 
feature due to our exphcit treatment of color is a color matrix element; the color factor for 
this diagram is 



' color 



1 ^ fA" A% 22 



33/2 



33/2 



(32) 



and is the same for diagram d2q and opposite for dig and ^25. The flavor and spin factors are 
specific to the reaction and will be discussed subsequently. The dig spatial overlap integral 
analogous to (A5) is 

1 



Ispaceidig) = 2\ d-'xd-'y -K(\x-y 



d'^ad^bd^c 0^(25 - A) 0^(26 - B) (pci^^c- C) 



5{A~d-B + h) {c\J{x)\a) {hc\J{y)\Q) , 

where the matrix elements of the currents (with F = 7°, 7*,!) are 

1 m„ 



(33) 



{q'\J{x)\q) = + 



(27r)3 



a g*(i?'-i?)-a; 



Uglg'TUgs 



(34) 
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{qVixM 



m„ 



,i{q'-q)-x 



(27r)3 



and 



{qq\Jiym 



1 m. 



i g«(9+9)-J/ 



Uqs^^qs 



(27r)3 

With these substitutions, and introducing a Fourier transformed kernel 



/C(g) = jd^xe'^-''K{r) 



we find 



(35) 



(36) 



(37) 



rf'^a d^c 0A(2a - A) 0^(2a - 2C - B) <pc{2c- C) 



Mb^^rvcse /C(|a-c| 



(38) 



and there are imphcit momentum constraints h = a — C and c = C* — c from the spatial 
integrations in (33) and the meson wavef unctions. 

As in the ^Pq model we take the nonrelativistic limit of this interaction for our decay 
amplitude, although it would be interesting in future work to investigate the effect of keeping 
the full relativistic amplitude (38). We also speciahze to the rest frame, ^4 = and B — —C. 
These substitutions give our final general result for the diq overlap integral, 

1 



' space 



(dlq) = 



(27r)= 



d^ad'^c (l)A(2d) 0^(2a + S) (f)c(2c + B) 



UhsJ^Vcso /C(|a-c| 



(39) 



The overlap integrals associated with the three remaining diagrams are given in Appendix B, 
and the detailed evaluation of a decay amplitude M.ls and a decay rate is illustrated in 
Appendix C for the decay p — > tttt with a j°Kj° interaction. 
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C. JKJ A— >^BC matrix elements: explicit results for SHO wavefunctions. 



The complete decay amplitudes for the channels spanned by our representative set of 
hght two-body meson decays are given below. These include decay amplitudes from 1) 
linear scalar, 2) constant scalar, 3) Coulomb OGE and 4) transverse OGE terms; a common 
factor is removed. 



Mls = 



(40) 



and the transverse OGE amplitude is shown as separate 4) 6ij and 5) —QiQj/Q^ contribu- 
tions. These five results for each channel are 



2^5 / b 



/ 1 3 A 4 / 1 5 A 

iFii-^;^;g + 77iFii-^;o;^j 



311/2 

T 
■33 

? 
■33 

23 



7T 



1/2 ( ^ 





X 



-|-— q;^ X 



iFi(^;|e)-^iFi(^;|e) 

/13 a 10 /I 5 A 
^^H2'2'^i + V^H2'2'V 

iFi(^;^;e)-^iFi(-;-;e) 



(41) 



^(3p,_.lSo + ^So) 



+ 



22 5^2 / b 



33 

210 



X 



313/251/2 
2^ 



71 



1/2 



So 



X 



3351/2 

2^ 



iFi(^;|e)-^iFi(i;|e 



15a 8 / 



2' 2 



8 /17 
45^^n2'2'^ 



■ 3351/2 



325V2 



as X 



17 /I 3 A 4 ^ /I 5 A 8 ^ /I 7 A 
iFi(2' + 3^^H2' 2'^i + 27^^n2' 2'^>' 

iFi(^;^;e)-^iFi(-;-;e)--iFi(-;-;e) 



; 2;C) 



2' 2 



(42) 



15 



23/25I/2 . 5 



35/2 \p2 

■3651/2 V/? 

23/2 



X 



Fi( - ^; + :^iFi( - ^; + ^iFi( - ^; 



1 5 



2' 2' 



15 



2' 2' 



225 



X 



a. X 



35/251/2'-^ 

25/2 



35/251/2 

23/2 
■33/251/2 



as X 



as X 



^2' 2 

•1 3 



^2' 2' 
1 5 



2' 2' 
1 7 



Fi (2 ; 2 ; - ^1^1 (9 ' 2 ' ^) " ^^^^(2 ' 2 ' ^) 



2' 2' 



27 



^2' 2' 



1 7 



2' 2' 



(43) 



_^(3P.-.3S. + ^So) 



2^5 / h 



(- 

37/2 \^^2 

29 



lFi( - |; - ^iFi( - I + 4i^i( " I 



H ^TT 

35 

25 



35/2 
26 



35/2^^ 

25 



33/2 



a. 



lFi( 
lFi(- 



9' 



lFi(-;^;t^;C 



1 .N 10 ^ / 11a 28 ^ / 1 3 A 
~2' ~ y - 2' 2' V + yi^il - 2' 2' 

~2' V + gi^il - 2' 2'^i ~ 2' 2' V 



(44) 



_^(3P,-.3S,+1S„) 



23/25 / 6 



/ 13 \ / ±o \ / ±( \ 



1 5 



1 7 



37/2 



2' 2' 



15 



2' 2' 



225 



2' 2' 



37 

23/2 
■3^ 



a, X 



16 



25/2 



37/2 
23/2 

■357^ 



a. X 



1 5 



^2' 2 



27 



^2' 2' 



27 



iFi(^;^;e) 



.~;.(3Po-iSo+iSo) 



29/^5 / b 

217/2 



39/2 
29/2 



iFi 



2' ~2' ~ y - 2' 2' + Ts'^n - 2' ^ 



32 

29/2 



32 
29/2 



1 

2'" 






2' 2'^' 


1 

2' " 






2' 2'^) 


1 

2'" 


"2'^) + 9^^^(~2'2'^) 


-^.F,( 
9 V 


2' 2'^' 



^(lp,_3s,+lSo) 



^ 29/^5 ( b 



37/2 ^^2 



Fi( - I - - I' h^) + lli^i 



2' 2 
2 



3 5 
2' 2 



35 

29/2 



-X 



35/2 
29/2 



"33/2"^ 

29/2 

a. 



33/2 



lFi(- 
lFi(- 



2'-y 9 ^- 2' 2'"^ ■ 27 

+ |iFi( - ^; ^; - ^iFi( - ^; |; C) 



2' 2 
1 

2' 2 
3 

2' 2 



1 1 .N 4 / 1 1 .\ 8 ^ / 13 
2;ej + 3iFi(--;-;ej-3iFi(--;-;^ 

^ A 16 / 11a 32 / 1 3 A 
"2' " T^^n - 2' 2' V + ^i^^l - 2' 2' 
L 1 
2"' 2^ 



22 5 / 6 



3^2 V/?2 



X 



iFi( - ^; + ^iFi( - I; + ^iFi( - h 
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D. JKJ decay rates: numerical results. 

To estimate the numerical importance of the scalar and OGE decay mechanisms we first 
specialize to p ^ tttt and consider each contribution in isolation. Using (6), (40) and (41), 
and the flavor factors in App.A, the decay rate due to linear scalar pair production alone 
would be 



-isKs 



TT 



-1/2/ 



1 3 



1 5 



2'2' 



45^ 



2'2' 



(49) 



The analogous decay rates due to Coulomb OGE and transverse OGE are 



iFi(|;|;e)-|iFi(^;|;e) 



n 2 



(50) 



and 



(51) 



These rates are shown in figure 5 for the parameter set as — 0.6, h — 0.18 GeV^ and 
rUq — 0.33 GeV. Evidently the dominant decay mechanism is the sKs interaction, pair 
production through the scalar confining potential. At j3 = 0.4 GeV the p tttt width 
predicted by the sKs model alone is 330 MeV, about twice the experimental 151 MeV. In 
comparison, transverse OGE gives a width of 3.9 MeV, and the color Coulomb interaction 
gives only 0.36 MeV. (The transverse and Coulomb OGE contributions add constructively, 
so the total width from OGE alone would be 6.7 MeV. These both interfere destructively 
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with the dominant sKs amphtude, so the total width we find from all three amplitudes is 
r^^^^ = 243 MeV.) 




Figure 5. The decay rate Tp^^^j^ assuming only sKs, j'^Kj"'" or j°Kj*^ decay interactions. 
Parameters a, = 0.6, 6 = 0.18 GeV^, = 0.33 GeV. 

Although the sKs interaction is usually found to be dominant in our representative set of 
decays, the OGE contributions are often comparable to sKs and cannot generally be ignored. 
To illustrate this, in figure 6 we show the numerical decay amplitudes we find in each decay. 
The normalization incorporates phase space, so these amplitudes squared give the physical 
decay rates. The magnitude of the experimental amplitude is indicated by horizontal lines. 
(The hi amplitudes are not shown because they are essentially identical to the bi amplitudes 
times a fiavor factor of y/S.) 
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Figure 6a-i. Decay amplitudes from sKs [br only), OGE (transverse) and OGE (coulomb) 
decay mechanisms, with parameters f3 = 0.4 GeV, b = 0.18 GeV^, = 0.6 and rriq = 
0.33 GeV. The decay amplitude from a constant scalar '^Pq term of 5*0 = —0.5 GeV is shown 
for comparison. 

The decay amplitudes due to Coulomb and transverse OGE are also shown in the his- 
tograms. The Coulomb gluon term is small in all the decays we have considered. This could 
have been anticipated; in the limit of a constant kernel K the Coulomb gluon decay ampli- 
tude actually vanishes, because the transition operator is then proportional to the fermion 
number operator Q = J d^xj^ (squared) and hence cannot pair produce. With a slowly 
varying kernel the Coulomb decay amplitude is nonzero, but remains small. Transverse 
gluon exchange in contrast is sufficiently large to make an important contribution in some 
channels. The most notable of these are Kq — ^Ktt and /o tttt, in which the OGE decay 
amplitude actually dominates the nonperturbative sKs amplitude (see figures 6h,i). Recall 
that the widths of scalar qq states are problematical in the ^Pq model because of a node in 
their decay amplitudes (see figure 2). The large additional OGE decay amplitude insures 
that the scalar qq states will be broad resonances, despite the node in the sKs or ^Pq decay 
amplitudes. 

Note that the overall scale of the total decay amplitudes predicted by the sKs+OGE 
decay mechanisms is too large relative to data by about a factor of two. It is nontrivial that 
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there is even approximate agreement; recall in contrast that the overall scale of the decay 
rates is not determined in the ^Pq model, but is fitted using the parameter 7. Since we regard 
the values as = 0.6, b = 0.18 GeV^ and ruq = 0.33 GeV as reasonably well established for 
light quarks, we doubt that the discrepancy reflects our choices for these parameters. For 
this reason we do not show an optimized fit to the experimental rates, which would give 
reduced values for ag or b. 

We suggest that the theoretical overestimate of scale may indicate the presence of impor- 
tant relativistic corrections to our nonrelativistic amplitudes, such asm/ E factors in (38); if 
included these would reduce the overall scale of rates and have little effect on the successful 
relative rate predictions. Other possibile sources of this discrepancy include our choice of 
qq wavefunctions (the overall scale of rates is very sensitive to the SHO width parameter /3) 
and our assumption that confinement may be treated as the exchange of a scalar quantum 
with a linear kernel in configuration space but no q-q-scalar form factor. 

A large negative constant 5*0 of up to Ri —1 GeV is often included in potential models, 
and is needed to give a best fit to spectroscopy. Of course it may not actually be present 
in the scalar potential S{r), and may instead simply be a way of subtracting off a fictitious 
constituent quark mass contribution of 2mg ~ 0.7 GeV from the qq spectrum. This constant 
would contribute a decay amplitude identical to a "^Pq coupling of 7 = (2'^/3^/^)5'o/mq; for 
completeness we have included such a constant term in our amplitude histograms, with an 
intermediate magnitude of Sq — —0.5 GeV. If an Sq of approximately this magnitude were 
actually present it clearly would make an important contribution, as it is larger than the br 
term and opposite in sign. We find however that any 5*0 of comparable magnitude leads to 
unrealistic D/S ratios in bi — > cutt and oi — > pir. In both decays D/S passes through zero 
near Sq = —0.3 GcV and diverges near Sq = —0.5 GeV with conventional parameter values 
of as = 0.6, p = 0.4 GeV, b = 0.18 GeV^ and rriq = 0.33 GeV, whereas S'o = leads to 
reasonable D/S values. Since the D/S ratios are inconsistent with a large negative constant 
5*0 we will assume Sq = subsequently. 
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E. D/S ratios in JKJ decay models 



In Sec.I.C we noted that the ratio of S- and D-wave amphtudes in the decays bi — > un and 
Oi — > pn allowed sensitive tests of the angular quantum numbers of the qq pair produced in a 
strong decay, and that the ^Pq model is in reasonable agreement with experiment, albeit for 
the rather large value oi P — 0.45 GeV. Recall also that the ratio of D/S ratios (ai/bi) in the 
^Po model is predicted to be —1/2 (neglecting minor phase space differences), independent 
of the Lqq = and Lgg = 1 radial wavef unctions. 

We find that the sKs interaction leads to very similar predictions to the ^Pq model. The 
sKs predictions for D/S ratios can be read from the decay amplitudes (41-48); for bi 
it is 



UJTT 
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(52) 



This rather complicated result is shown as a function of (3 in figure 7, together with the 
^Po prediction (21) and the corresponding oi — > fm results. The sKs and ^Pq decay models 
evidently give remarkably similar D/S ratios, and are essentially indistinguishable in these 
reactions. 

The ratio of ai/6i D/S ratios with the sKs interaction is predicted to be 
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(53) 



just as in the ^Pq model. Thus, to the extent that the ^Pq model is successful in explaining 
D/S ratios, this is a success of the sKs model as well. Some similarity could have been 
anticipated since both the ^Pq and sKs decay models assume pair production with ^Pq 
quantum numbers, and differ only in the presence or absence of spatial correlations with 
incoming quark lines. 
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Figure 7. D/S ratios in bi — >■ cuvr and Oi — >■ /ovr predicted by the '^Pq model (dashed) and the 
sKs model with S{r) = br (solid). 

The simple ratio in (53) follows in the j^Kj" model as well, and appears to be generally 
true if the initial lines have no spin-flip amplitude and identical (ai,6i) and (p, ci;,7r) radial 
wavefunctions are assumed. In contrast, the j"'"Kj"'" transverse OGE interaction has an 
initial-line spin-flip amplitude, so that the ratio departs from —1/2. This may be useful as 
a signature of the OGE component in the physical decay amplitude. 

Although these D/S ratios are often cited as an argument against an OGE decay mecha- 
nism, the actual OGE calculation does not appear to have been carried out in the literature. 
(We emphasize that OGE is not equivalent to the "'^Si" decay model, for which results do 
exist P,p!0|.) We find that the full (Coulomb + transverse) OGE D/S ratios are given by 
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These are shown in figure 8; evidently the hypothesis of OGE dominance of these decays 
can indeed be rejected, because the OGE decay mechanism predicts D/S amphtude ratios 
of the wrong sign for both Oi pir and bi — > un. 




Figure 8. D/S ratios in bi ujtt and ai —>■ pir assuming pure OGE decay amphtudes. 

In reahty both sKs and OGE decay amphtudes are present, and figure 6 shows that 
the OGE contributions are not neghgible here. It is especially interesting to investigate the 
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combined effect of sKs and OGE amplitudes, because it may be possible to identify the 
individual contributions through interference. In the following we set M(ai) = M{bi) = 
1.23 GeV and M{p) = M{u!) = 0.78 GeV to remove small phase space effects. First in 
figure 9 we show the D/S ratio in bi ujtt as a function of the OGE coupling strength as 
for a realistic sKs amplitude (b = 0.18 GeV^), for several values of the wavefunction scale (3. 




Figure 9. D/S ratio in bi un with combined sKs and OGE decay amplitudes. 

Evidently the effect of OGE is to reduce D/S somewhat, and for the conventional as = 0.6 
this leads to a more realistic value of the wavefunction scale, (3 ~ 0.35 GeV. Recall that a 
rather large jS ~ 0.45 GeV was required to fit the observed bi and ai D/S ratios with the 
sKs interaction treated in isolation. 

A more definitive test of the presence of OGE contributions involves the ratio of ai/bi 
D/S ratios; this is exactly —1/2 (assuming identical phase space and spatial wavefunctions) 
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for any sKs or ^Pq decay strength. We show this ratio versus the OGE couphng in figure 10, 
for the same parameter set as figure 9. 




Figure 10. The ratio of D/S ratios, (ai —>■ p7r)/(6i t^Ti"), with combined sKs and OGE 
decay amphtudes. 

For as = 0.6 we expect a reduction in the ai/&i D/S ratio to ~ —0.46 to —0.47. The 
current experimental ratio of —0.35(9) suggests a larger departure from the sKs and '^Pq 
prediction of —1/2, but is consistent with either theoretical result within errors. A more 
accurate experimental determination of these D/S ratios would be interesting as a test of 
the expected OGE corrections to nonperturbative sKs or '^Pq decay amplitudes. 
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III. SUMMARY, CONCLUSIONS AND FUTURE APPLICATIONS. 



In this paper we have discussed the current theoretical understanding of open flavor 
strong decays and proposed new "microscopic" quark and gluon mechanisms for these de- 
cays. We began by reviewing the ^Pq decay model and introduced a diagrammatic repre- 
sentation for its decay amphtudes. We then carried out detailed calculations of the decay 
amplitudes for a representative set of light meson decays, first in the ^Pq model and then 
in microscopic decay models, which assume that the qq pair is created by an OGE mech- 
anism or by nonperturbative pair production from the scalar confining potential. These 
results were derived analytically using well established light quark interactions and SHO 
wavef unctions. 

We find that with conventional quark model parameters the light qq open-fiavor meson 
decays are usually dominated by qq pair production from the scalar confining potential. We 
refer to this as the sKs (scalar-kernel-scalar) decay model. This explains the success of the 
phenomenological ^Pq decay model, since the 'ipip scalar current produces pairs in a ^Pq state, 
which is assumed a priori in the ^Pq model. The dimensionless pair production strength 
7 is not determined theoretically in the ^Pq model, and is treated as a free parameter. In 
contrast, in our calculations the absolute decay rates are completely determined by the qq 
wavefunction scale /3 and known QCD interaction parameters as, b (string tension) and rrig. 
In the sKs decay model the ^Pq strength 7 corresponds approximately to the dimensionless 
combination h/mq(5. 

Although there are differences in detail between the ^Pq and sKs models due to the 
different overlap integrals, we find that the numerical predictions for relative decay rates 
and amplitudes in the decays we have considered are remarkably similar. In addition to 
explaining the success of the ^Pq model, our description of decays also accounts naturally 
for the absence of "hairpin" diagrams; since the microscopic pair production interactions 
produce qq pairs in color-octet states, hairpin diagrams are forbidden without additional 
interactions. In the usual ^Po-like decay models one assumes pair production of a qq color 
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singlet, and hairpin diagrams are simply ignored without theoretical justification. Other 
tests of the color state of qq pair production in decays would be useful to discriminate 
between the ^Pq model and our color-octet description. 

An unsolved problem is that the overall scale of the amphtudes in the full sKs+OGE 
microscopic decay model is too large by about a factor of two, given conventional quark 
model parameters. This discrepancy of scale may be due to our neglect of relativistic effects 
such as the niq/E external line normalizations in (38). 

Although we have shown that OGE qq pair production is usually dominated by the non- 
perturbative sKs decay amplitudes, wc do find numerically important OGE contributions 
in some channels, notably ^Pq — > ^Sq+^Sq. In that case OGE pair production actually 
dominates the nonperturbative sKs amplitude, and insures broad /o(1300) and Kg(1430) qq 
resonances. OGE pair production also makes characteristic contributions to some observ- 
ables such as the D/S amplitude ratios in bi — > cutt and oi — > fm. It may be possible to 
identify OGE contributions in these and similar multiamplitude decays. 

In future work it would be interesting to apply these microscopic decay calculations 
to charmonium (where the transverse OGE contribution should be much smaller) and to 
strong decays of light baryons, for which new data is expected from the CEBAF experimental 
program. The simple exercise of extending these light meson calculations to higher-L^g initial 
mesons would also be very useful, as it appears likely that large departures from ^Pq and 
sKs model predictions due to OGE contributions could be found among the many decay 
channels available to these states. 
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APPENDIX A: DIAGRAMMATIC FORMULATION OF THE ^Pg MODEL. 

1. General results. 

One may simplify calculations of decay amplitudes in the ^Pq model by developing a 
diagrammatic description. This is nonessential for '^Pq model calculations of meson decays 
since there are only two diagrams, but in our generalized decay models there are four dia- 
grams, and it is useful to distinguish their contributions to prove relations between them. 
We anticipate that the diagrammatic description will also be useful in the more complicated 
combinatorics of baryon decays. 

We begin by noting that the pair production component of the ^Pq Hamiltonian (2) can 
be written in terms of creation operators as 



We associate both the coupling constant g and the external-line spinor bilinear ^ ["^A-s'^-fo] 
with an effective ^Pq qq pair production vertex. There is no additional factor of — i, unlike 
conventional field theoretic Feynman rules, because we are determining the matrix element 
of Hi instead of the T-matrix. 

We assume nonrelativistic qq wavefunctions for the initial and final mesons; in our no- 
tation a meson state is of the form 



with implicit spin and flavor wavefunctions that are the usual nonrelativistic quark model 
forms. (For details of our conventions for the wavefunctions see reference ]15|], Appendix 
A.; note that a factor of l/(27r)'^/^ was inadvertantly omitted from the normalization of the 
wavefunction in equation (A15) in that reference.) 

One may now evaluate the Hamiltonian matrix element for the decay A— >BC in terms 
of straightforward overlap integrals. Schematically, we have a matrix element of the type 




(Al) 




(A2) 



{hhcc\ h\ d) 

^ ' ks ~K! 



aa) 



(A3) 
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and it is useful to distinguish two Feynman diagrams (figure Al) which we refer to as di if 
the produced quark goes into meson B and d2 if it goes into C. 





Figure Al. qq meson decay diagrams in the ^Pq decay model. 

These are drawn so the ordered quark and antiquark content of the state can be read from 
the diagram by drawing a vertical line through it. Thus, in di the initial state is |A) = \ad) 
but immediately after the pair creation we have 6j(i||aa) = \qqaa). This is then rearranged 
to \qaaq), which is projected onto the final state \BC) = \bbcc). The odd number of crossed 
lines in each diagram implies a (—1) phase (the "signature" of the matrix element) due to 
permutation of quark and antiquark operators. 

Specializing to diagram di, the associated matrix element is of the form 

{B C \ H [\ A"^ [l-^ -^signature ' -^flavor ' Ispin+space • (-^^) 

The spatial overlap integral associated with diagram di (before the spin matrix element is 
taken, which gives Ispin+space (c^i)) is 

Ispace(c?i) ^ ////// d^O: d^o. d^b d% d^cd^c 
(j){a - a) 6{A ~ a - a) (f)*{b -b) 6{B -b -b) (p*{c-c)6{C-c-c) 



meson A meson B meson C 

/ 777 r 1 — * — * — * — 

d'kg-^ u^^v_-,- 6{k-b)6{-k-^) 5{a-c)S0-b) (A5) 



^ V ' spectator lines 

pair production amplitude 

This result can be read directly from the diagram (see figure A2), using the pair production 
vertex (Al) and the fact that each unscattered "spectator" line gives a factor of 6{ki — kj). 
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Figure A2. Determination of the spatial overlap integral for diagram di 



Thus the decay amplitude is a 21-dimensional integral involving 21 one-dimensional delta 
functions. Only 18 of the integrations can be carried out trivially, which leaves a three- 
dimensional overlap integral times a 6{A — B — C) momentum- conserving delta function, 



^space space 



5{A-B-C 



as in (4). These overlap integrals are explicitly (setting A = and B = —C) 



spaceidl 



d'-^k (j)A{2k - 2B) (t)l{2k - B) (f)},{2k - B) ■ 

Ek 



d''k (t)A{2k + 2B) (f)l{2k + B) 0^(2A; + B) ■ g-^ 

Ek 



(A6) 

(AT) 
(A8) 



^ space ( ^2 ) 

The spin factor and labels Sq, Sq in these overlaps depend on the reaction being considered, 
and are determined by the external line labels attached to the diagrams, as illustrated below. 



2. An illustrative ^Pq decay: p —>■ vrvr. 

Here and in Appendix C we will use the decay p+(+i) tt+tt" to illustrate our tech- 
niques, and simply quote results for other cases in the text. For this decay the flavor states 
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are \p~^) = \n^) = —\ud) and |7r°) = {\uu) — \dd))/\^. Figure A3 shows the evaluation of 
the flavor factor for diagram di. 




Figure A3. Determination of the flavor factor of +1/^/2 for diagram di in p+ n^n°. 

The flavor factors associated with the other decays we consider can be determined sim- 
ilarly, and are tabulated below. For a given decay, such as /2 tttt, we give the flavor 
weight factors for each diagram in a speciflc charge channel (in this case /2 tt+tt^). The 
resulting decay rate for /2 tt+tt^ must then be multiplied by a multiplicity factor JF to 
sum over all flnal flavor states (here J-' = 3/2 for tc^tt^ and tt°tt°). 



Table Al. Flavor Weight Factors. 


Generic Decay 


Subprocess 


I flavor (di) 


Iflavor{d2) 




P — ^ TTTT 


p+ ^ 7r+7r° 


+1/V2 
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/ ^ TTTT 


/ — > TT + TT^ 


-1/V2 


-1/v^ 


3/2 


a — ^ pvr 


a+ p+7r° 


+1/V2 


-1/v^ 
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b ^ un 


6+ — s> u;7r+ 


+1/V2 




1 


h ^ pn 


h p+TT^ 


-1/V2 




3 


K* Kvr 


K*+ ^ K+7r° 


+1/V2 
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The "flavorless" decay amplitudes quoted in the text, equations (7-15) and (40-48), corre- 
spond to unit flavor factors, Ifiavor{di) = +1 and Ifiavor{d2) = ±1. (The sign is chosen 
so the contributions of di and d2 add rather than cancel.) To convert these to physical 
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decay amplitudes one should multiply (7) or (40) by Ifiavor{di) from Table Al, and multiply 
the resulting total decay rate by J^. The K* is a special case because only one diagram 
contributes; in this case (7) or (40) should instead be multiplied by Ifiavor{di)/2. 

The spin states are \p{+z)) = | t T ) and Itt) = (| t I ) ~ I i T ))/V2; taking the spin 
matrix element analogously to figure A3, we find a spin factor of [(— 1/2)^^ J for this 
diagram. Combining these results and using the explicit Dirac spinor matrix element in 
Appendix B, we find 

hfi{d\) — -^signature' ^flavor ■-^spin+space(^l) (-^^) 
(-1) (+1/V2) 

where 

W.pace(c?i) = —9 J d^kcPA{2k-2B)ct>*j,{2k-B) ct>*c{2k -B) . (AlO) 

The second diagram ^2 has an opposite fiavor factor, so the full result is 

— * — * 

Since the overlap integrals satisfy 7(^2, —B) — +I{di, B), hfi is odd under Qb — —^b, and 
can be written as 

hfi \/2 space (c^l) 

= +^ — [ d^k (t)A{2k - 2B ) (f)*B{2k - B ) 0^(2^ -B){k^ + iky) . (A12) 

y2 tTlq J 

We have also substituted ruq for Ef. in (AlO) to recover the nonrelativistic limit, which gives 
the ^Po model. 

Frequently in the literature these decay amplitudes are evaluated using SHO qq wave- 
functions, which leads to closed-form results. For Lgg — mesons our Gaussian momentum 

space wavefunction (p^kq — kg), with k = kq = —kg, is 

<i>m = e-'y^^\ (A13) 
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We note in passing that the Lgg — 1 wavefunction used in this paper is 

o3/2 1 h. 

'^-(2^^) = ^ I e-^V2.^ . (AM) 

— * 

On substitution into (A12), (f){2k) gives a n'^n" decay amphtude of 

and an A4io ^'^^^ amphtude of 



Substitution into (6) then gives the total decay rate 

which with the identification 7 = g/'^mq is the usual ^Pq result. 

Many features of this decay rate could be anticipated on general grounds. Specifically: 

1. The {g/mqY dependence follows from the 0{g) of the amplitude and the nonrelativistic 
limit of (2). 

2. The factor of Mp comes from phase space and the relation EbEc/Ma — Mp/A. 

3. The threshold dependence is expected for a P-wave final state. 

4. Dimensional arguments then require a factor of 1//?^. 

5. An exponential in [P/ (3Y is expected given SHO wavef unctions. 

Only the overall numerical coefficient and the factor of 1/6 in the exponential (present in 
all 3Po SHO A— >BC qq meson decays) require detailed calculation. 
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APPENDIX B: OVERLAP INTEGRALS AND SPIN FACTORS. 



In the text we discussed the evaluation of only one of the four decay diagrams present 
in A— >BC meson decay in JKJ models. Here we give the overlap integrals associated with 
all four diagrams: 



d^ad^c (j)A{2d) 0B(2a + B) 0c(2c+ B) 



{2n) 



Ubsb^Vcs, /C(a-C) UcsJUas, 



(Bl) 



(with implicit momentum constraints specific to this diagram of 6 = a + B and c — —B — c), 

Ispaceidlg) = J J d\d% M^^) 0b(26 " B) (f>*c{^d + B) 



(with c = —B — c, d — —a and b — B — b), 

1 



space{d'2q) — ~l~ 



(27r)= 



^•^0 d^b (l)A(2d) 0^(26 - B) (/)^(2a - B) 



(B2) 



K,{d-b) UbsJ^u 

aSa 



(with c — d — B and b — B — b), 



'space(c?2,-) = J J d'ad'c 0^(2a) cj>*B{2d-B) cj>*c{2c+ B) 



(B3) 



ic{d — c ) 



(B4) 



(with a — —d, b = B — d and c = —B — c). 

In our evaluation of decay matrix elements in JKJ models we also require spin matrix 
elements, which involve the nonrelativistic 0{p/m) matrix elements of Dirac bilinears with 
r = 7°, 7 and I and Pauli spin matrix elements for various spin states. These are 
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lim 



2mq 



lim 

■y/c— *0 



r = 70 

{q + q' )5ss' - i{s'\a \s) x {q' - q) T = 7 

r = i 



_L.Jss\a\0) ■ {q + q) T = 7° 
(ss|a |0) r = 7 

^{ss\a\0)-iQ-Q) r = I 



(B5) 



(B6) 



lim 

v/c-^O 



VgsTVgfs' 



5ss' r = 7° 

{q + q')Sss' -i{s'\a\s) X {q' -q) T = ^ 

r = i 



1 

2to, 



(B7) 



The explicit scattering and pair-production matrix elements of Pauli spinors in terms of 
the spherical basis vectors e± = ^{x + iy)/\/2 and Cq = z, are 



(T 1^1 T) = -(i 1^1 i)=eo 
(T |a| i) = V2e_ 

(i \a\ T) = -V2e+ 

(Tl^lT) = -(II<?II) = -eo 
(T|a|I) = x/2e+ 
(I|a|T) = -x/2e_ 
(TT|^|0) = V2e_ 
(Tl|5^|0) = (iT|^ |0) = -eo 
(iI|5'|0) = V2e+ 



(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 
(B15) 
(B16) 
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APPENDIX C: p ^ TTvr IN JKJ DECAY MODELS 

We will illustrate a JKJ decay calculation in detail using the reaction p'^{+z) tt+tt", 
as in our discussion of the ^Pq model in Appendix A. We will evaluate the p — > tttt decay rate 
using the j'^Kj° interaction in (27-29), and simply quote results in the text for the other cases. 
For each type of JKJ interaction we find that the four diagrams give equal contributions, so 
we discuss only diagram dig, shown below, in detail. 




B 



Figure CI. Decay diagram dig. 



iq 



For all the JKJ interactions we consider this diagram has a signature of (—1), a color 
factor of +2^/3^/^, and a flavor factor (for p+ vr+vr") of +1/2^/^. Differences in the 
relative importance of the decay mechanisms arise in the spin-space overlap integrals and 
the r spinor matrix elements. 

In the text we derived the spatial overlap integral for diagram dig in the nonrelativistic 
limit, equation (39); 

1 



<■ space 



(dlq) 



{2ny 



d^ad^c (pAC^a) (/)^(2a + S) (f)c{2c + B) 



UbsJ^Vcs-, lC{\d-c\) 



^CScJ^^aSa 



(CI) 



To evaluate this for p^{+z) tt+tt", note that for the color Coulomb interaction F = 7° 
the spinor bilinears give 
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2mn 



{sbSc\a\0) ■ (b + c) 5s,sa 



(C2) 



Attaching the spin wavefunctions to the diagram diq gives an overall factor of (—1/2) and 

sets Sb =1 and Sg = I, so the spin matrix element is {(j)gg = (— 1/2)(| ||a|0) = — e+/\/2. 
With the substitutions Q — {a — c) — {b + c) and E = (a + c)/2, hfi becomes 



hfi{d^,)^ (-!)( + 





1 X 


) (+33/2)( 


+ 21/2 J 



signature color 



flavor 



II 



(27r)3 2mg 



ic{Q) 0p(2E + g) 0;+(2E + g + s) 0;„(2s-Q + 5) (C3) 



Without additional information about the spatial wavefunctions this is the final result. To 
proceed we assume standard quark model Gaussian wavefunctions (A13) for the mesons. 
The d^T, integral can then be carried out, which gives 
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/ d'Q {a)g,-QIC{Q) exp { - ^(Q - (C4) 



where we have again abbreviated P/ (5 — x. 

We may evaluate this intermediate integral / by returning to coordinate space using (37) 
and integrating over d^Q, which gives 



I = 33/V/2/33 / {a)gg ■ {]-B + -f3'f ) K{r) e i^'^" 



(C5) 



For {a)qq in this reaction we have 



/27r 



(C6) 



and the angular integrals can then be carried out using a spherical harmonic expansion of 
the plane wave in (C4). This gives 

3 



I ^ -2^/^371^(3^ YuinB) / r^drK(r)e- 

Jo 



(C7) 
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Since the kernel is simply K{r) — ag/r in this case, the radial integrals can be evaluated in 
terms of confluent hypergeometric functions using the general formula 



f 

Jo 



where = {n + i + l)/2. This result suffices for all the overlap integrals we encounter in 
JKJ decay models with SHO wavefunctions and power-law kernels. 
Applied to the Coulomb case, our result for the integral / is 

I = -2^/V«./?^ xIfA |;^) - |,F,(i; l;^)] e-^'"'Y,,{ns) , (C9) 



where ^ — Using this integral we obtain our final result for diagram dig, 



= "^^-'"^^'"^ - |iFi(^; e-^''^-Y,,{ns) . (CIO) 



On multiplying by 4 for the four equivalent diagrams, we obtain the full result for hfi for 
p^{-\-z) — > 7r+7r°. We can abstract the j'^Kj^ decay amplitude M. ls in (6) from this; removing 

the p"*" — >• tt+tt" flavor factor of l/\/2, we have 



3d 



iFi(^;|;e)-|iFi(|;|;e) 



-0^2/12 



(Cll) 



which is equivalent to the result quoted in (41). Using (6), we find that the total p — >■ tttt 
decay rate due to color Coulomb OGE (treated in isolation) is 



X 



iFi(^;|;e)-^iFi(^;^;e) 



1 5 



(C12) 



Alternatively, for the specific case of a Coulomb kernel we can evaluate the momentum 
space integral (C3) directly. The required integrals and others needed to evaluate decay 
amplitudes for P-wave quarkonia are given below. 

r3/2 r .13 



/ d^Q (l/g^) e-«(Mo)^ = ^ 2,F,(i,|,agl 



(CIS) 



(C14) 
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7^0 



3/2 f 1 9 ^7 1- 

-^-k iFi{^,^,aQl)+^QoiQoj iFi(-,-,aQ2)| e-»Q§ , (C15) 



aV2 i 3a 



^3/2 /- 2 3 7 -» 

^Yy^-jy^ {SijQok + SjkQoi + SkiQoj) iFi(-, -, aQo) 

+^aQoigo,-Qofe iFi(Hag^)| g-^^o . (C17) 
This approach apphed to p — > tttt leads to the form 

which is equivalent to (C12). 

Relations between confluent hypergeometric functions with different indices allow these 
results to be written in various forms. We typically express our final results for decay ampli- 
tudes as linear combinations of confluent hypergeometric functions with constant coefficients 
and a common first index a, times a centrifical factor of x^^'^ . These rearrangements are 
straightforward using the recurrence relations 

iFi (a; c; x) = l^^^y i^i (a - 1; c - 1; x) + j^^iFi (a - 1; c; x) (C19) 

and 

XiFi(a;c;x) = ■^^-^^^^^^y^|iFi(a - 1; c - 2; x) - iFi(a - 1; c - 1; x) | . (C20) 

The decay rates due to transverse OGE and the confining interaction may be derived 
using the same techniques, with the minor complication that transverse OGE also has initial- 
line spin-flip contributions. These results are quoted in Section II. D. 
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